Preprint Oxford: 
OUTP-93-30P 
Preprint IFUNAM: 
FT94-39 Feb/94 



On 
ON. 

Effect of wavefunction renormalisation in N-flavour 
£ '■ QED3 at finite temperature 

M • by 

I.J.R. AitchisonQ and M. Klein-KreislerQ 

Department of Physics 
. Theoretical Physics 

1 Keble Road 
Oxford 

Os 1 0X1 3NP 



Submitted to: Physical Review D 



Abstract 



a,: 
a: 

> 
X 

A recent study of dynamical chiral symmetry breaking in N-flavour QED3 at finite tem- 
perature is extended to include the effect of fermion wavefunction renormalisation in the 
Schwinger-Dyson equations. The simple "zero-frequency" truncation previously used is 
found to lead to unphysical results, especially as T — > 0. A modified set of equations is 
proposed, whose solutions behave in a way which is qualitatively similar to the T = solu- 
tions of Pennington et al. [5-8] who have made extensive studies of the effect of wavefunction 
renormalisation in this context, and who concluded that there was no critical N c (at T=0) 
above which chiral symmetry was restored. In contrast, we find that our modified equations 
predict a critical N c at T 7^ 0, and an N — T phase diagram very similar to the earlier 
study neglecting wavefunction renormalisation. The reason for the difference is traced to the 
different infrared behaviour of the vacuum polarisation at T = and at T 7^ 0. 
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1 Introduction 



In a recent paper with collaborators [1], we studied dynamical chiral symmetry breaking in 
N— flavour QED 3 at finite temperature, in the large N approximation. Using an approximate 
treatment of the Schwinger-Dyson equation for the fermion self-energy, we found that chiral 
symmetry was restored above a certain critical temperature, which itself depended on N. 
The ratio r of twice the zero-momentum and zero-temperature fermion mass to the critical 
temperature turned out to have a value of about 10 (approximately independent of N), 
which is considerably larger than a typical BCS value, but consistent with previous work [2] 
using a momentum-independent self-energy. We found evidence for a temperature-dependent 
critical number of flavours, above which chiral symmetry was restored, and the N — T phase 
diagram for spontaneous mass generation in the theory was presented. A fuller account of 
some of the relevant details, and of the possible relevance to high T c superconductivity, is 
contained in [3]. 

The question of the existence, or not, of a critical N in analogous calculations at zero 
temperature is still to some extent controversial. The original calculations of Appelquist 
and co-workers [4] found chiral symmetry breaking only for N < N c where N c = 32/7T 2 , 
but this work made a possibly crucial appeal to perturbation theory (in l/N) to justify 
the neglect of wavefunction renormalisation and the use of a simple bare vertex. This step 
has been strongly criticized by Pennington et al. [5,6,7,8], who in a series of papers, using 
Schwinger-Dyson equations with increasingly elaborate non-perturbative vertices satisfying 
the Ward and Ward-Takahashi identities, have found no evidence for any N c - rather, the 
fermion mass simply decreases exponentially with N. On the other hand, an alternative 
non-perturbative study by Atkinson et al. [9] suggested that chiral symmetry is restored at 
large enough N; but this paper is also criticized in recent work of Curtis et al. [8]. Finally, 
Kondo and Nakatani [10] have examined the effect of imposing an infra-red cut-off on the 
Schwinger-Dyson equations, including wavefunction renormalisation and various ansatzes 
for the vertex. For the Pennington- Webb [5] vertex, which we also shall use, Kondo and 
Nakatani [10] obtained a (cut-off dependent) N c for the case of a "large" infra-red cut-off, 
and N c — > oo for the case of a "small" cut-off. 

Our previously mentioned calculations at finite temperature [1] made precisely the same 
perturbative approximation as Appelquist et al. [4], by neglecting wavefunction renormalisa- 
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tion and using the bare vertex. The arguments and results of Pennington et al. [5-8] certainly 
provide strong motivation to go beyond that approximation, and investigate whether the con- 
clusions of our finite temperature study survive a better treatment of the vertex. The present 
paper reports the results of such an investigation, using a simplified (but not perturbative) 
form of the vertex, but otherwise following [1] as closely as possible. In brief, we find that 
the conclusions of [1] are, in fact, essentially unchanged, though the formulation of a simple 
extension of the formalism of [1] to include wavefunction renormalisation turns out to be 
not completely straightforward. It seems that the crucial ingredient in obtaining a finite N c 
is the softening of the infra-red behaviour at finite temperature. Our calculations therefore 
suggest a natural physical interpretation of the "large" infra-red cut-off regime introduced 
phenomenologically (at T = 0) by Kondo and Nakatani [10]. 

2 Simple Schwinger-Dyson equations at T ^ 

We shall choose a simplified ansatz for the vertex, and we begin by introducing that choice 
within the context of the zero-temperature S-D equations; then we shall pass to the finite 
temperature case. In Euclidean space, the S-D equation for the fermion propagator is 

S~ F \p) = Sf-\p) - e 2 / ^t*M*)A/«,(?)HM (1) 

where the superscript denotes the bare quantity, Sp 1 (p) = (1 + A(p))$ + Y,(p) is the inverse 
fermion propagator, and q = k — p. We shall not make any change here in our previous (and 
common) choice for A^ u - namely, it is approximated by the sum of massless fermion bubbles. 
We shall also continue to work in Landau gauge. Turning then to T u , the Ward-Takahashi 
identity may be written as 

q u T»(k,p) = S F 1 (k)-S F 1 (p). (2) 
Taking the limit q u — > yields the Ward identity 

^ = r (P , P ) (3) 

OPu 

which ensures that the full vertex is free of kinematic singularities. Both (2) and (3) should 
hold to all orders of perturbation theory. 
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Ball and Chiu [11] have given a vertex which satisfies both these relations, but has an 
unconstrained transverse part (which, however, is believed to be unimportant in the infrared 
region [11, 5-7]). In fact, a more general vertex including not only the full Ball-Chiu vertex 
but also a non-trivial transverse part has recently been studied by Curtis et al. [8], at zero 
temperature. This vertex satisfies (2) and (3) and also correctly reproduces the leading 
asymptotic behaviour known from perturbation theory. For our purposes, however, the 
salient fact is that the results (as to dynamical mass generation, N c etc.) of this most recent 
work are in qualitative agreement with earlier studies [5,6] where a much simpler choice of 
vertex was made (albeit one no longer satisfying (2) and (3)), namely 

r»(k, P ) = Y(i + A(k)). (4) 

Similar approximations in four dimensions have been shown to give results which are in 
fairly good agreement with those obtained using the full vertex. It seems reasonable to hope 
that such agreement will persist in our finite temperature case, and we shall therefore now 
adopt the Pennington- Webb [5] vertex (4). 

Inserting (4) into (1) and separating the scalar and spinor parts we find (see also Eqns. 
(2.35) and (2.36) of [3]) 

M(v) = -^- 1 [d*k 1 M{k) (5) 

{P) 4ttW 1 + A(p) J q 2 + U(q) k 2 + M\k) K> 

and 

f*, iP 2 -k 2 f-{q 2 f 1 



A ^ 167T 3 7Vp2 J dk q 2 [q 2 + Tl(q)} k 2 + M 2 (k) ^ 
where a = e 2 N is understood to be fixed as N varies, II is the vacuum polarization, and we 
have introduced the mass function 

M ^ - ttwy (7) 

The explicit appearance of the 1/N factor in (6) is, of course, the reason for the claim 
that to leading order in 1/N wavefunction renormalisation can be ignored and A set to 
zero. But Pennington and co-workers have convincingly argued [5-8] that, at least at the low 
momenta relevant to dynamical mass generation, the integral can provide a compensating 
iV-dependence (as after all happens for £ or M), so that A is by no means of order 1/N. 
These authors [5] have obtained numerical solutions to equations equivalent to (5) and (6), 
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which show no sign of a critical N c . We follow them in retaining (6), and pass now to the 
finite temperature case, following the usual prescriptions. As in our previous work [1] (see 
also [3]) we define 

p = (P , p), P = |p|, Po = (2m + 1)tt/(3 ((3 = l/k B T) 

k = (ifeo, k), K = \k\, k = (2n + l)ir/(3 (8) 

Q = (go, q), Q=|q| = |k-p|, q = 2(n-m)7r/(3. 

Integrals over the temporal component of a fermion loop momentum are replaced by 
infinite sums over odd Matsubara frequencies, while bosonic loops are evaluated by summing 
over even frequences. The vacuum polarisation and the fermion functions A and £ become 
functions of (the modulus of) the momentum, and of the temperature, and acquire a discrete 
index n corresponding to the Matsubara temporal component. Since we are here concerned 
with the effect of introducing the non-perturbative vertex (4), we shall follow [1] and [3] in 
retaining only the fj, — v — component of the photon propagator A^ u , and ignore all but 
the zero-frequency (n = 0) component, so that 

A^(q ,Q,f3) = 5^o 5„ /[Q 2 + MQ,P)l (9) 
where, to an excellent approximation [1], 



H.W./J) = § 



\Un2 ( -Ti \ 

QP + _ exp (— QP) 



(10) 



In the same spirit, we shall also ignore the frequency dependence of the kinematic factors 
in (6). In this zero-frequency limit, then, both A and S become independent of the Matsubara 
frequency index, and the sums over these indices in (5) and (6) can be performed explicitly 
so as to yield the following simple equations for the temperature-dependent mass function 
M(P,(3) and wavefunction renormalisation A(P, (3): 

M(PR , * 1 fj2 V M(K,0) tanhfy^ + .M^) 

{ ,P ' 8Nn* 1 + A(P, (3) J -Q^ + U (Q,P) ^K* + M*(K,(3) {> 

and 



With regard to equations (11) and (12), we note that the first (for M) is the same as 
that for E in ref . [1] , except for the coefficient in front of the integral which acquires a factor 
(1 + A(P ) f3))~ 1 in the present case. In [1], of course, A was set to zero. We also note that 



A does not appear under the integral in (12), but is given by a simple integral involving Ai. 
This latter feature is a consequence of the ansatz (4); in general we would have had to deal 
with two coupled integral equations for M. and A. 

It would seem that all that now remains is to solve equations (11) and (12) numerically. 
However, it turns out that there is an unsuspected problem (or so we regard it) with Eq.(12), 
as we now explain. 



3 Problem with the sign of A in the solution of (11) 
and (12). 



It is a standard result in zero-temperature field theory (see for example [12]) that the com- 
plete wavefunction renormalisation factor Z = 1 + A satisfies the relation 



Indeed, the angular integral in the zero-temperature expression (6) for A was evaluated 
analytically in [5] and found to be negative definite (which clearly implies A < 0), and the 
corresponding numerically evaluated Z satisfied (13). Unfortunately, the same results do not 
hold for our approximate expression (12) for A at finite temperature, as we shall now see. 

Consider the angular integral in (12), namely 



where Q 2 = P 2 + K 2 — 2PK cos 0, and we have included the (dimensionless) factor Kj a 
coming from the two-dimensional phase space d 2 k. As it is not possible to evaluate / an- 
alytically, we have had to resort to numerical evaluation. Clearly it is difficult to give a 
complete picture of / as a function of all three variables P, K and f3, but it turns out that 
for wide ranges of these variables I is predominantly positive, although it is negative in the 
region around K m P. These features are illustrated in Figure 1, which shows / versus K/at 
for two values of P (one "large" compared to the natural scale a, the other small), and two 
values of (inverse) temperature (3. It is clear from (15) that / is always negative at the point 



< Z < 1, 



(13) 



which implies 



- 1< A < 0. 



(14) 




(15) 
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K = P, and that the quantity IIo acts as a kind of regulator for the 1/Q 2 singularity asso- 
ciated with the photon propagator. In particular, from Fig. 1 we see that the width of the 
region where / is negative decreases with increasing (3 (decreasing T). These latter details 
are understandable from the form of il , Eqn.(lO), which has a temperature-independent 
linear term (aQ/8) and another which depends on temperature as (aC/8(3)exp(—0Q/C) 
where C = (16&i2)/tt. 

It is of course difficult to read off from Fig.l what the eventual sign of A in (12) will 
be, since much might depend on the relative weighting attached to the region K P, 
and to the remainder, in the if-integration of (12). Nevertheless, it does seem likely that, 
especially at the low temperatures characteristic of dynamical mass generation, the relative 
insignificance of the negative parts of / in Fig.l will imply that A in (12) turns out to be 
positive, contradicting (14). This is indeed the case. 

We have solved (11) and (12) by an iterative procedure, as follows. In zeroth order, we 
took = and — E, the solution of Eq.8 of our previous paper [1] (which is the same 
as (11) with A = 0). Inserting into the right hand side of (12) gave the first iterate 

A^\ This was then substituted into the coefficient in front of the integral in (11), and the 
latter evaluated using as the input function, so as to yield the first iterate M^. This 
was substituted back into (12) to give A^ 2 \ and so on. The procedure was continued until 
convergence to within 2 % was achieved. Note that, as in [1], we work with a momentum 
cut-off at A = a, and scale all momenta, temperatures and masses by a. 

For the sake of illustration we consider the case N = 1. We have obtained stable and 
converged solutions to Eqns (11) and (12), with properties we now describe. First, the 
behaviour of M(P, (3) is qualitatively similar to that of E(P, f3) found in [1] - namely, M. 
is constant for P/a < 1CT 2 , and falls rapidly to zero for larger P/a. Further, the zero- 
momentum value of Ai starts to fall rapidly when (3a goes from 2000 to 1000, indicating 
the possibility of a finite critical temperature. Indeed a plot of Ai(0,T)/a versus ksT/a 
suggests a critical temperature of order ksT c ~ 10~ 3 a:, quite similar to the T c 's found in [1]. 
However, all these results involve an A which is greater than zero, and hence a Z violating 
(13). Fig. 2 shows the corresponding solution of (12) for N = 1 as a function of momentum, 
for different temperatures. The behaviour of the full Z — 1 + A is quite different from 
that found by Pennington and Webb [5] from the zero-temperature equation (6) - which, as 
mentioned above, always satisfies (14). 
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We might wonder whether for very low temperatures our solution to (12) goes negative, 
but this does not happen. Instead, as the temperature decreases the value of A(0,/3) rises, 
approaching a value of approximately 0.6. This behaviour bears no resemblance to the 
zero-temperature result of Pennington and Webb [5] . 

It is certainly possible that the condition (14), which relies on unitarity in Minkowski 
space, may not necessarily hold in the Euclidean space appropriate to T ^ 0. There may, for 
example, be "heat-bath" creation processes for temperatures above the pair-creation thresh- 
old which would cause a violation of (14). But we shall take the view here that the regime 
relevant to dynamical mass generation is definitely a low temperature one (ksT <C a), and 
that consequently we should hope to find an A which is negative and qualitatively similar 
to that of [5]. After all, the essential aim of the present work is to examine the effect of a 
reasonable temperature-dependent extension of [5] on the existence or otherwise of a critical 
N c . We therefore reject the solutions of (11) and (12) described above, and seek a modified 
Eqn (12) which will give a temperature-dependent A satisfying (14), for the (low) tempera- 
ture with which we are concerned. 



4 Modification of the "A" equation to secure 

-1<A<0 at Tyo. 

We want to understand why the zero-temperature expression (6) for A satisfies (14), while 
our approximate finite-temperature version of it gives A > 0. We believe the answer lies in 
our dropping all but the zero frequency components in II, E, and A. First, it is clear that if 
all components are kept, the finite temperature equations must correctly reduce to the zero- 
temperature ones as T — > 0. More particularly, by retaining only the zero component we 
have effectively lost a dimension (compare (12) with (6)), and this is crucial for the following 
reason. Suppose we consider the zero temperature limit of (15), in which n — > aQ/8. In 
this case the integrand of (15) is dominated by very large positive values associated with the 
regions close to = and 2tt, whereas for intermediate values of the integrand is negative 
but very much smaller. This is why, as noted earlier, / is predominantly positive. However, 
if we were simply to multiply the integrand of (15) by sin0, so as to mimic the three- 
dimensional phase-space at zero temperature, we would effectively eliminate the unwanted 
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large positive contribution, and enhance the negative ones. In fact, we have checked that 
introducing such a factor by hand in (15), and integrating from to n, renders I negative 
for all P and K - and hence ensures A < 0. Hence we believe that a proper "reconstruction" 
of the full phase space, at least near T = 0, would succeed in changing the sign of A as 
desired. 

Unfortunately, it seems a formidable task to attack the fully coupled equations, including 
all frequency components. Instead, we shall seek here a simple modification of (15), which 
behaves qualitatively in as similar a way as possible to the zero-temperature kernel in (6), 
but is temperature-dependent. In this way we hope at least to model the effect which the 
Pennington- Webb vertex (4) would have in a more realistic frequency-dependent calculation. 

The modification we propose is motivated in part by the notion that the approximation 
of retaining only the zero frequency components is best justified at high temperatures. In 
this (small (3) limit, the quantity n in (10) reduces to the (temperature-dependent) value 
(2a£n2/n(3) - but of course we are really interested in low temperatures. However, dynamical 
mass generation is a low-momentum phenomenon, and furthermore the integrands in (5) or 
(6) are certainly enhanced for q 0, so that we may hope that it may be reasonable to try 
a "small j3Q" approximation, rather than simply a "small f3" one. In the small (3Q limit, n 
of (10) reduces again to (2a£n2/ir(3), plus 0((3Q) 2 corrections. These considerations lead us 
to explore the result of replacing n in (12) by a constant, A 2 say, which we will take to be 
an adjustable parameter with a value of order a 2 or less. We will choose it so as to obtain 
a wavefunction renormalisation as much like ref. [5] as possible. 

In replacing n in (12) by A 2 we are of course also altering the Q — > behaviour of 
the kernel, making it less singular. To "compensate" for this we might think of dropping 
the K in the phase-space factor, at the same time as replacing n by A 2 . To explore these 
possibilities we shall study (a) the "exact" angular integral (with the K factor included) of 
(15), called I(P,K,f3); (b) the modified kernel with n — > A 2 and no X-factor, namely 



and (c) the kernel with n — > A 2 but retaining the X-factor, i.e. the quantity K x 7a- 
Of these we shall prefer the one which, on the one hand, best captures the low-momentum 
behaviour of I(P,K,f3) (since this is the relevant region for dynamical mass generation), 
and on the other gives a wavefunction which most resembles the qualitative behaviour of the 




(16) 
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T = results of ref. [5]. 

It turns out that I a can be integrated analytically, with the result 



_-2tt 



P 2 - K 2 1 [(P 2 - ^ 2 ) + A 2 ] [(P 2 - K 2 ) - A 2 ] 



1- J L + 



(17) 



A 2 A 2 {[(P - K) 2 + A 2 ][(P + K) 2 + A 2 ]p_ ' 
Figure 3 shows a plot of (17) as a function of P and (scaled by a), for the case A 2 = a 2 . 
We note that, most importantly, Ia is almost always negative (in contrast to / of (15)) except 
for a small region of large K and P values, which is not the important region for dynamical 
mass generation. In fact, the behaviour of I a is similar to that of an analogous integral which 
arose in an approximate treatment [6] of the three-dimensional zero-temperature equation 
(6): namely, a rapid variation with K for K < P, and a much slower variation for K > P, 
the sign being (in [6]) negative throughout. These features are true of I a for a substantial 
range of A 2 , and ensure A < as we shall see in the next section. 

Figure 4 shows a comparison of J, I a and K x I a for P/a — 6.851 x 10~ 4 and A 2 = a 2 . It 
is remarkable how well the low K behaviour of / is captured by I a, and how poorly K x I a 
performs. The large positive contribution in / for larger values of K is, as noted above, an 
undesirable feature and the reason for obtaining A > in section 3. On the other hand, if 
too small a value of A 2 is chosen (for example A 2 of order 0.1a; 2 ) we find that I a becomes 
too large and negative, with the consequence that the corresponding A begins to approach 
-1. From the appearance of the factor (1 + A)" 1 in (11) we expect that such A's will be 
associated with too large a value of M, contradicting the basic assumption M <C a. The 
above considerations already strongly suggest that Ia will be the most satisfactory kernel at 
T 7^ 0, and we now proceed to discuss the results of using this kernel in place of I in (12); 
we shall also briefly describe the results of using K x Ia- 

5 Results using a modified "A" equation 

We have solved the original Eqn (11), together with one or the other of the following modified 
equations for A: 



MP,P) = t£t5 dKI A (P,K) = l J (18) 



16^7r 2 io " ' ^K 2 + M 2 (K : /3) 



a r tanh ! \l K 2 + M 2 (K, (3) 

A K ^P^) = - KW - 2 dKKI A (P,K) ; V ==== , (19) 

16Nn 2 Jo ^/K 2 + M 2 (K,(3) 
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with I a given by Eqn. (17). We followed the same iterative procedure as described in 
Section 3, obtaining stable and converged solutions for A and E. The results using Aka 
can be quickly summarised. Figure 5 shows A K a(P,P) versus P/a for A 2 = a 2 and N = 1; 
we note that Aka is negative, as required, but very small in magnitude. Clearly this is a 
consequence of the very small (and negative) value of K x I a as seen in Fig. 4. Figure 6 
shows the corresponding A4(0,T)/a versus ksT/a, for N — 1. This figure is very similar to 
the A = 0, N = 1 case shown in fig. 3 of ref. [1], as would be expected from the small value 
of A K a shown in Fig. 5. In particular, the ratio r = 2M(0, 0)//cbT c remains close to 10. 
Qualitatively, the effect of the small negative Ika is to produce a small upward shift of the 
Ai(0,T)/a curve, relative to the A = case; this is easily understood as being associated 
with the (1 + A)^ 1 factor in (11). We have obtained the corresponding phase diagram using 
Ika, but it differs only rather slightly, and predictably, from fig. 7 of ref. [1]. 

The results using I a are more interesting. For A 2 ~ 0.1a 2 , we found that A A was positive 
for a small range of P, near P w a; otherwise, A A (P) was negative as required. More 
seriously, for this value of A 2 we found that M.(P) was only about one order of magnitude 
smaller than a, so that the condition M. <C a characteristic of dynamical mass generation in 
QED 3 (and of [5]) was no longer true (implying that the approximation of II by the massless 
bubble would need re-examination). For A 2 = a 2 , however, we did obtain values for M. and 
Aa which were in much closer qualitative agreement with the zero-temperature results, and 
we now concentrate on the results for this value of A 2 . 

In Figure 7 we show Ai(P, (3) as a function of momentum, for N = 1 and different /3's, 
while Fig. 8 shows the corresponding A A 's, all satisfying (14). These A A 's are all very 
reasonable-looking, reminiscent of those at T = [5,6]. The rapid decrease of M as (3 de- 
creases from (3 = 500/q; to 130/a suggests the existence of a critical temperature near 10 _2 o;. 
In Fig. 9 we show Ai(0,T) as a function of temperature, from which it appears that the 
critical temperature T c is such that ksT c « 7.8xl0~ 3 a. Figures 7 and 9 can be compared 
with Figures 2 and 3 of [1]; the conclusion is that the effect of including wavefunction renor- 
malisation in the way described is to increase both the mass M. and the critical temperature 
T c by roughly a factor of two; again, this can be qualitatively understood as a consequence 
of the (1 + A)^ 1 factor in (11). Indeed, the ratio r (= 2A^(0,0)/fc B T c ) comes out as 10.27, 
almost identical to the value for iV = 1 which we quoted in [1] without using wavefunction 
renormalisation. 
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Fig. 8 shows that the dependence of A A on (3 is not very strong, disappearing altogether 
at the high momentum end. In Fig. 10 we show A A at P = 0, as a function of T (note that 
(18) actually has non trivial solutions even for M = 0, so that we can investigate the region 
T >T C for A A ). It can be shown easily that A A (P = 0, T) -> as T -> oo. 

Finally, we come to the main question we originally set out to answer - the existence of 
a critical iV c , or otherwise. Fig. 11 shows AA(0, (5) versus N, for various fixed temperatures. 
As in our previous work [1] , we are not able to obtain reliable results for M. /a much below 
1CT 5 . Nonetheless, as before, it seems reasonable to infer that M. does vanish beyond some 
finite N c which itself depends on T, increasing as T decreases. We can extrapolate the 
Ai/a curves to find the critical values N C (T) for the various temperatures, and thus obtain 
the phase diagram shown in Figure 12, and which is very similar to the one we obtained 
previously with A = [1]. 

Comparing Figure 11 with Figure 5 of [1], we notice that for a given temperature, the 
critical number of flavours has increased in the case where A ^ 0: for example, at (5 — 10 4 /a, 
N c was just above 1.8 in the A = case, whereas it now has a value slightly above 2.2. 
Although probably too much should not be made of such a relatively small difference, it is in 
fact simple to understand it, following an argument of Pennington and Webb [5]. Comparing 
our (11) above with (8) of [1], we can see that a measure of the effect of A ^ can be obtained 
by replacing N in the A = case by N(l + (A)) where (A) is an average value for A. Hence 
a solution of the A = equations for a given Na=o will effectively correspond to a solution 
of the A ^ equations with the identification Na^o(1 + (A)) = Na=o- Since (A) < 0, the 
critical number of flavours for a given T increases in going to the A ^ case. 

Fig. 13 shows A A as a function of P for various values of N, at fixed T. Although we 
have not been able to explore as wide a range of A's as Pennington and Walsh [7] , this figure 
is qualitatively similar to their Figure 3, for N ~ 2, which encourages us to think that our 
modified ^A-equation, (18), is physically reasonable. However, it is interesting to note that 
in the zero-temperature results of [5] and [7], A(P = 0) seems to approach -1 with increasing 
N, suggesting that the corresponding N e ff = (N A=0 )/(1 + (A)) becomes so large that any 
sign of criticality in iV (which might have been present in the A = case) disappears. Even 
for smaller TV's, where (A) is considerably different from -1, the results of [5]- [8] do not show 
any indication of a sudden decrease in the dynamical mass at some critical N, such as is 
seen in our Fig. 11. 
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Nevertheless, it should be stressed that we are not able to obtain reliable numerical 
results for temperatures below ~ lCT 4 a, so that we are not able to say what the precise 
zero-temperature limit of our N C (T) might be - if indeed it exists at all. In addition, we 
must not forget the inexact nature of the "A 2 modification". In short, it is quite possible 
that N C (T) — > oo as T — > 0, which would be in agreement with the conclusion of Pennington 
et al. [5-8]. 

If this suggestion is correct, we need to understand why our non-zero temperature results 
still point so clearly to a critical N, even though we have included a non-zero A. The 
answer seems surely to lie in the important alteration which the finite-temperature vacuum 
polarisation Ho makes to the infrared regime of the fermion self-energy. As we have discussed 
earlier, the Q — > limit of Ho in Eqn. (10) is a finite temperature-dependent term 2a£n2/7i[3 
- a phenomenon called "thermal screening". The "A 2 modification" we used in equation (18) 
for Aa also captures this softening, for small Qf3. By contrast, the zero-temperature IT of 
Pennington et al. [5-8] (inherited from Pisarski [13]) behaves as aQ for small Q and has 
no infrared softening. We believe that it is this infrared screening, associated with the 
temperature-dependent n , which is crucial for the existence of an N C (T), just as it was 
found to be vital to the generation of a large value for r [1]. 

6 Conclusions 

The work of Pennington et al. [5-8] strongly indicated that the existence [4] of a critical num- 
ber of flavours N c , above which chiral symmetry was unbroken in QED3 at zero-temperature, 
was an artefact of incorrectly ignoring wavefunction renormalisation (via an unjustified ap- 
peal to perturbation theory in l/N). We have been interested in extending the study of chiral 
symmetry breaking in QED 3 to finite temperature. In [1] we found clear evidence for the ex- 
istence of an N c , dependent on T, but we made the approximation of ignoring wavefunction 
renormalisation. The present study has been aimed at removing that approximation, and 
studying the effect of including wavefunction renormalisation on chiral symmetry breaking 
in QED 3 at T ^ 0. 

The simplest generalisation of the model of [1] to include wavefunction renormalisation, 
in which we adopted the "zero-frequency" approximation, turned out to lead to unphysical 
results, in that they failed to show any similarity to the zero-temperature results of [5-8] as 
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T — > 0. We regarded this as reason for discarding that model, and adopted instead a modified 
equation for the wavefunction renormalisation at T 7^ 0, which gave results consistent with 
the T = case. We found, using the modified model (namely the kernel I a of sections 4 and 
5), that although wavefunction renormalisation was now included in a way qualitatively very 
similar to that of [5] - [8], nevertheless at finite temperature we still found clear evidence for 
an N c , in contrast to the results of [5]- [8]. The essential reason for the difference seemed to 
be the characteristically different infrared behaviour of the vacuum polarisation, at zero and 
at finite temperature. Indeed, the inclusion of wavefunction renormalisation, in the manner 
described in Section 4 above, gave results very little different, overall, from those of [1]. In 
particular, a value of about 10 for the ratio r seems to be remarkably robust, while the values 
of .M(0, 0) and T c depend more on the model for A. 

The conclusion that it is the infra-red behaviour which is crucial in obtaining a finite N c is 
supported by the calculations of Kondo and Nakatani [10]. Using (among other vertices) the 
Pennington- Webb vertex of equation (4), these authors found that with an infra-red cut-off e 
of order 10 _4 a, a fermion mass was generated for N < N c , where N c depended logarithmically 
on e. Very roughly, we might interpret e as corresponding, in our calculation, to the "thermal 
mass" (n ) 1//2 , which is of order T 1 / 2 . We then have a physical interpretation of the perhaps 
somewhat artificial cut-off introduced by Kondo and Nakatani [10]. 

These authors [10] also found that the behaviour of the dynamically generated mass for 
N near N c , using the vertex of equation (4), was consistent with mean field theory: that is, 
Ai ~ (N c — iV) 1 / 2 . By contrast, in our previous paper [1] we presented evidence (see figure 
6 of [1]) which suggested a behaviour of the form M exp[— C/(N C — TV) 1 / 2 ]. However, 
we have re-examined those calculations and are now less confident that a firm conclusion 
can be drawn regarding the behaviour near N c , where the mass is very small and numerical 
accuracy increasingly difficult to control. In fact, a more conservative interpretation of both 
our present and previous data is that they are quite consistent with mean field theory, at 
least in the regime (for N not too close to N c ) where the calculations are most reliable. 

The present calculation can, of course, be criticised for invoking the somewhat ad hoc 
"A 2 modification" for the A equation (though we stress again that this modification does 
the job required of it, namely to produce sensible looking solutions at low T). It would 
clearly be much more satisfactory to attack the equations for the coupled n^0 components 
of S, using the full n ^ components of the photon propagator as given in [3]. In this 
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way the true effects of the n ^ components, especially at low T, could be identified. No 
doubt a start on this problem could be made by going back to the simplifying approximation 
S(P,/3) « £(0,/3), which is likely to be quite reliable (see Eqn (2.54) of [3]). We hope to 
return to this problem elsewhere. 
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Figure Captions 

Fig. 1 Angular integral I(P,K,/3) at P = 6.85 x 10 _4 a, 0.809a and (3 = 1000/a, 100/a. 

Fig. 2 The wavefunction renormalisation as a function of scaled momenta, at various (scaled) 
temperatures, for N = 1. 

Fig. 3 Angular integral I A (P, A) for A 2 = 0.3. 

Fig. 4 Comparison of J, I A and K x I A for P/a = 6.851 x 10~ 4 and A 2 = a 2 . 
Fig. 5 A KA (P,(3) versus P/a for A 2 = a 2 and iV = 1. 

Fig. 6 The scaled self-energy m{T)/a (= Ai(0,T)/a) as a function of scaled temperature for 
N = 1 and A 2 = a 2 using A^a- 

Fig. 7 The scaled self-energy as a function of scaled momentum for N — 1 and (3a = 10 4 , 
2000, 500 and 130. The parameter A 2 is fixed at the value A 2 = a 2 . 

Fig. 8 A A (P, (3) as a function of scaled momentum for fixed N — 1 and (3a = 10 4 , 2000, 500 
and 130, with A 2 = a 2 . 

Fig. 9 The scaled self-energy m(T)/a (= M(0,T)/a) as a function of scaled temperature for 
N = 1 and A 2 = a 2 using A A . 

Fig. 10 A A (0,T) as a function of scaled temperature for iV = 1 and A 2 = a 2 . 

Fig. 11 The scaled self-energy m((3)/a (= M(0, (3)/a) as a function of iV at various scaled 
temperatures. 

Fig. 12 The phase diagram of QED 3 , in our approximation, with non-trivial wavefunction 
renormalisation. The critical line separates the region where there is dynamical mass 
generation (Ai ^ 0) from that in which there is not (Ai = 0). 

Fig. 13 The wavefunction renormalisation A A (P, (3) as a function of scaled momenta, for 
N = 1, 1.5, 2, 2.2 at pa = 10 4 . 
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